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The influence of a topological defect in graphene on the ground state of elec- 
tronic quasiparticle excitations is studied in the framework of the long-wavelength 
continuum model originating in the tight-binding approximation for the nearest 
neighbour interaction in the graphitic lattice. A topological defect that rolls up 
a graphitic sheet into a nanocone is represented by a pointlike pseudomagnetic 
vortex with a flux which is related to the deficit angle of the cone. The method 
of self-adjoint extensions is employed to define the set of physically acceptable 
boundary conditions at the apex of the nanocone. The electronic system on a 
graphitic nanocone is found to acquire the ground state condensate and current of 
special type, and we determine the dependence of these quantities on the deficit 
angle of the nanocone, continuous parameter of the boundary condition at the 
apex, and the distance from the apex. 
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1 Introduction 

Topological phenomena are of great interest and importance because of their uni- 
versal nature connected with general properties of the space. Topological defects 
in the quasirelativistic fermionic matter can induce vacuum quantum numbers. 
A general theory of the vacuum polarization by a pointlike topological defect of 
the vortex type in two dimensional quantum systems of massless Dirac fermions 



was elaborated in Refs. [H [2]. In the present paper we apply this theory to the 
study of the ground state polarization in graphene with a topological defect (see 
also Refs. [HI [4]). 

Carbon atoms in graphene compose a planar honeycomb lattice with one va- 
lence electron per each site. The primitive cell is rhombic and contains two atoms, 
thus the graphene lattice consists of two rhombic sublattices. The first Brillouin 
zone is a regular hexagon with corners corresponding to the Fermi points; among 
six of them, the two oppositely located ones are inequivalent. Electronic quasi- 
particle excitations in graphene are characterized by a linear and isotropic dis- 
persion relation between the energy and the momentum in the vicinity of the 
Fermi points, where the valence and conduction bands touch each other. Us- 
ing the tight-binding approximation for the nearest neighbour interaction in the 
honeycomb lattice, an effective long-wavelength description of electronic states in 
graphene can be written in terms of a continuum model which is based on the 
Dirac-Weyl equation for masless electrons in 2 + 1-dimensional space-time with 
the role of speed of light c played by Fermi velocity v ~ c/300 [3 El [7]. The 
one-particle Hamiltonian operator of the model takes form 



where ajL and a 2 ^ are the 4x4 matrices belonging to a reducible representation 
composed as a direct sum of two inequivalent irreducible representations of the 
Clifford algebra in 2 + 1-dimensional space-time. The one-particle wave func- 
tion possesses 4 components, which reflects the existence of 2 sublattices and 2 
inequivalent Fermi points (valleys). 

Unlike the conventional case of spinor electrodynamics in 2 + 1-dimensional 
space-time (see, e.g., Ref. [8]), the parity transformation in the continuum model 
of graphene implies the inversion of both spatial axes and the exchange of both 
sublattices and valleys |9], 



H (o) 



ihv (afojdi + a 2 {0) d 2 ) , 



(1) 




(2) 



where 



PH^ = -H^P, P 2 = i 
The time reversal implies the exchange of valleys [10J, 



(3) 




(4) 



where 



T(H®y = # (0) T, T 2 = -/. 



(5) 



The matrix of the spatial inversion can be presented as 



P = 2Si?, 



(6) 
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where ^ 

is the pseudospin, and R satisfies commutation relations 

[R,al 0) ]_=[R,a 2 (0) ]_ = Q (8) 

and exchanges the sublattice indices, as well as the valley indices. 

In the second quantization picture, one can consider ground state expectation 
values: 

the P-condensate 
and the i?-current 



p(x) = (vac|^ f (x)P^(x)|vac) (9) 



j(x) = (vac\^(x)cxR^(x)\vac), (10) 

where x = (vt, x 1 , x 2 ), a = (a 1 , a 2 ), and |vac) denotes the ground state (vacuum). 
Evidently, quantities (9) and (10) are vanishing in the case of Hamiltonian given 
by (1), which corresponds to the idealized strictly planar graphene with 
all interactions neglected. In reality, the layers of graphene are corrugated at 
mesoscopic scales [TTJ [T2J fTH] . and namely the effects of curvature in graphene 
samples are addressed in the present paper. Therefore, our starting point is the 
ground state expectation value of the time-ordered product of fermion fields in 
the form 

(vac|T^(x)^(y)|vac) = (x\(hv^V '^{y), (11) 

where \I/ = ^7°, and V M (/i = 0, 1, 2) is the covariant derivative in curved 2 + 1- 
dimensional space-time. Restricting ourselves to static backgrounds (Vo = <9 ) 
and using Eq.(ll), we get 

p(x) = -itr(x\P(ihvd - Hy x \x) (12) 

and 

j(x) = — itr(x\atR(ihvd — i/) _1 |x), (13) 

where 

H = -ihva ■ V (14) 
is the Dirac-Weyl Hamiltonian on a curved surface, 

[at, aJ'] + = 2gV'l, (15) 

and Qjji is the metric of this surface. Further, using the Wick rotation of the time 
axis, Eqs.(12) and (13) are recast into the form which exhibits explicitly their 
time independence, 

p(x) = --tr(x|Psgn(if)|x) (16) 
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and 



j(x) = -~tr(x\aRsgn(H)\x). 



(17) 



In the present paper we compute the P-condensate and the P-current in 
graphene with a topological defect. 

2 Topological defects 

Topological defects in graphene are disclinations in the honeycomb lattice, re- 
sulting from the substitution of a hexagon by, say, a pentagon or a heptagon; 
such a disclination rolls up the graphitic sheet into a cone. More generally, a 
hexagon is substituted by a polygon with 6 — Nj sides, where Nj is an integer 
which is smaller than 6. Polygons with > (A^ < 0) induce locally positive 
(negative) curvature, whereas the graphitic sheet is flat away from the defect, as 
is the conical surface away from the apex. In the case of nanocones with > 0, 
the value of Na is related to apex angle S, 



and Nd counts the number of sectors of the value of 7r/3 which are removed from 
the graphitic sheet. If N d < 0, then —N d counts the number of such sectors which 
are inserted into the graphitic sheet. Certainly, polygonal defects with N d > 1 and 
N d < — 1 are mathematical abstractions, as are cones with a pointlike apex. In 
reality, the defects are smoothed, and N d > counts the number of the pentagonal 
defects which are tightly clustered producing a conical shape; graphitic nanocones 
with the apex angles 5 = 112.9°, 83.6°, 60.0°, 38.9°, 19.2°, which correspond to 
the values N d = 1, 2, 3, 4, 5, were observed experimentally [H]. Theory predicts 
also an infinite series of the saddle-like nanocones with — A^ counting the number 
of the heptagonal defects which are clustered in their central regions. Saddle- 
like nanocones serve as an element which is necessary for joining parts of carbon 
nanotubes of differing radii and for creating Schwarzite [15] , a structure appearing 
in many forms of carbon nanofoam [IB]. As it was shown by using molecular- 
dynamics simulations [T7j, in the case of N d < —4, a surface with a polygonal 
defect is more stable than a similarly shaped surface containing a multiple number 
of heptagons; a screw dislocation can be presented as the N d —>■ —oo limit of a 
6 — Ad-gonal defect. 

The metric of a conical surface with a pointlike apex has the form 



where r and <p are polar coordinates centred at the apex, and — oo < rj < 1. The 
intrinsic curvature of the cone possesses a 5 2 (x)-singularity at its apex, vanishing 
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at x 7^ 0, and parameter 77 enters the coefficient before this singularity term. 
Quantity 2irrj for < rj < 1 is the deficit angle measuring the magnitude of the 
removed sector, and quantity —2tti] for —00 < i] < is the proficit angle mea- 
suring the magnitude of the inserted sector. In the case of graphitic nanocones, 
parameter 77 takes discrete values: 



rj = N d /Q. 
Using Eqs.(15) and (18), one gets 

a r = a\ 0) , = (1 - ^V" 1 ^, 

and the Dirac-Weyl Hamiltonian on the cone takes form 

H = -ihv {a( )<9 r + «^ ) r 1 [(* ~ v)' 1 ^ - * s ] ] 
The second-quantized fermion field operator is presented as 



(19) 



(20) 



(21) 



E 



dE\E\ 

h 2 v 2 



exp [-iEx°(hv) x ] ip En (x.)a En + 



+ 





^ / ^^[-^AH- 1 ]^^, (22) 



where Z is the set of integer numbers, a En and ag n (b En and b En ) are the fermion 
(antifermion) creation and destruction operators satisfying anticommutation re- 
lations 



a En, a Eh 



S(E-E) 



I — " nn 5 

tJ + v EE 

and V'Bn(x) is the solution to the stationary Dirac-Weyl equation 

Hi/j En (x) = Ei/; En (x). 
The ground state is defined conventionally by relationship 

OE n |vac) = &e n |vac) = 0. 



(23) 



(24) 



(25) 



Solutions to the Dirac-Weyl equation form a complete set and are orthonormal- 
ized in a way which is usual for the case of the continuum 



2n 



J dip J drVsVinW^EnW = 2^ 



5{E-E) 

u 1 — u nni 

V EE 



(26) 



where g = detgjji = (1 — r]) 2 r 2 , and a factor of 2 in the right hand side of the last 
relation is due to the existence of two inequivalent Fermi points (valleys). 

As it was shown in Ref. f3], the fermion field on a graphitic nanocone obeys 
the Mobius-strip-type condition: 

if?(vt, r, Lp + 2vr) = - exp(-i3irr]R)^(v t, r, ip), (27) 

where r\ is given by Eq.(19). Condition (27) in the case of odd Nd involves the 
exchange of sublattices, as well as valleys. Note that since R 2 = I, the exchange 
is eliminated after double rotation 

W(vt, r, ip + 47r) = cos(N d 7r)^(vt, r, cp); (28) 

that is why the mention of the Mobius strip seems to be relevant. 

By performing a singular gauge transformation (see Ref. [3] for more details), 
one gets the fermion field obeying usual condition 

r, ip + 2vr) = -W{vt, r, <p), (29) 

in the meantime, Hamiltonian (21) is transformed to 



H' = —ihv I cnjL d r + o3 \r 1 



(30) 



Thus, a topological defect in graphene is represented by a pseudomagnetic vortex 
with flux A^7r/2 through the apex of a cone with deficit angle Note that, 

due to commutation relations 

[P, R}- = [P, E]_ = [T, R]+ = [T, S] + = 0, (31) 

discrete symmetries of spatial inversion and time reversal are maintained: 

PH' = -H'P, T(H'Y = H'T. (32) 



3 Solution to the Dirac— Weyl equation 

Vacuum expectation values are independent of the matrix representation used, 
therefore a choice of representation is a matter of convenience. As it was already 
noted, the ahy and a^-matrices (and, consequently, S) are of the block-diagonal 
form. Since the i?-matrix satisfies relation (8), it can be unitarily transformed to 
the block-diagonal form also: 

URU^ 1 — ( q — I ) ' ^«(o)^ _1 = a(o). (33) 



6 



Thus Hamiltonian attains the block-diagonal form after this unitary transforma- 
tion: 

H" = UH'U- 1 =( H q h ! ) ■ (34) 

To be more precise, let us assign the definite sublattice and valley indices to 
components of the initial fermion field in the following way [18J: 

^ = (V A+ , ^b+, ^a-, #b-) T , (35) 

where subscripts A and B correspond to two sublattices and subscripts + and — 
correspond to two valleys. After performing the singular gauge transformation 
and the unitary one, we get \J r " with components mixing up different sublattices 
and valleys. The appropriate solution to the Dirac-Weyl equation takes form 

i>En = (^En,U ^En-lf \ (36) 

where the two-component functions satisfy equations 

H s ^En,s = E^ ErltS , S = ±1. (37) 

Corresponding to Eq.(35), the ahy and a^-matrices can be chosen in the 
form 



, (J 2 \ 2 / a 1 , 

"(o) =- I o a 2 j , «(o) = I o -a 1 ' ' , ,N 



where a 1 and a 2 are the off-diagonal Pauli matrices. Then the matrices of spatial 
inversion and time reversal in the initial representation take form 

p=fv;y T=i(y n ). ( 39) 



v ° / ' V 1 , 

Separating the radial and angular variables in the solution to Eq.(37), 

gmAr) ei{n+i) , ), (40) 



we get that the radial components satisfy equations 

Dl >s \ ( f En>s (r) \ _ f f En , s (r) , M | , 

D n , s ) I g En Jr) r ' ' 



where 



D n , s = hv[-d r + r 1 (1-T]) 1 (sn-r])], 

Dl s = hv [d r + r -1 (l -r 1 )- 1 {sn+ 1 -2^)] . (42) 



Let us consider graphitic nanocones with 1 > 77 > — | and 77 = — 1, and define 
quantity 

p = I [I - 1^(77) + V] (1 - I?) -1 , 1 > »7 > -5 fa 7^ 0), (43) 

A pair of linearly independent solutions to Eq.(41) is written in terms of the 
cylinder functions. In the case of 1 > 77 > | (Nd = 5, 4, 3), the condition of 
regularity at the origin is equivalent to the condition of square integrability at 
this point, and this selects a physically reasonable solution. Thus, in view of the 
orthonormality condition (26), the complete set is given by regular modes with 
sn > 



fEn,s( r ) \ = 1 ( Ji(i-r,)-i- F (kr) 

9En,s{r) ) 2yV(l -77) V sgn(E) J i(1 _„)-i +1 _ F (A;r) 

and regular modes with sn < 

fEn,s{r) \ = 1 ( Jj/(i_,)-i +F (At) 

#En, s (r) y 2 v /tt(1 — 77) V -sgn(^) Ji'(i_„)-i_i+F(^) 



/ = sn, (44) 



I' = -sn, (45) 



where k = \E\(hv) 1 , and J M (w) is the Bessel function of order /z; note that F is 
integer belonging to range 5 > F > 1 in this case. 

In the case of \ > n > (N d = 2, 1), > 77 > -\ (N d = -1, -2, -3) and 
77 = — 1 (iV d = —6), there is a mode, for which the condition of regularity at 
the origin is not equivalent to the condition of square integrability at this point: 
both linearly independent solutions for this mode are at once irregular and square 
integrable at the origin. To be more precise, let us define in this case 

^ = I § [sgn(77) - 1] , \ > n > -\ (j) + 0), (4g) 

Then the complete set of solutions to Eq.(41) is chosen in the following form: 
regular modes with sn > sn c 



( fEn,s{r) ^ _ 

V 9En,s{r) J 2^(1-77) 
regular modes with sn < sn c 

fEn,s( r ) \ = 1 

9En,s{r) ) 2^(1-77) 



Ji(i-n)-i-F(kr) 
sgn(E)Ji {1 _ v) -i +1 _ F (kr) 



Ji>0--T,)- 1 +F(kr) 
-sgn(S) J,/(i_„)-i_i +f (At) 



, / = s(n — n c 



(47) 



, l' = s(n c -n), 
(48) 
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and an irregular mode 

fEn c ,s( r ) \ _ L 



X 



9En c A r ) J 2y/n(l - rj) [1 + sin(2^) cos(Ftt)] 

x f sm(v E )J- F (kr) + cos(is E )J F (kr) \ , . 

y sgn(£) [sm(i/ B ) Ji- F (kr) - cos(v E ) J-i +F (kr)} J ' "' 

note that F belongs to range < F < 1 in this case. Thus, the requirement of 
regularity for all modes is in contradiction with the requirement of completeness 
for these modes. The problem is to find a condition allowing for irregular at 
r — ► behaviour of the mode with n = n c , i.e. to fix v F in Eq.(49). To solve 
this problem, first of all one has to recall the result of Ref . [19] , stating that for 
the partial Dirac Hamiltonian to be essentially self-adjoint, it is necessary and 
sufficient that a non-square-integrable (at r — > 0) solution exist. Since such a 
solution does not exist in the case of n = n c , the appropriate partial Hamiltonian 
is not essentially self-adjoint. The Weyl-von Neumann theory of self-adjoint op- 
erators (see, e.g., Ref. |20j) is to be employed in order to consider a possibility 
of the self-adjoint extension for this operator. It can be shown (see Ref. [3]) that 
the self-adjoint extension exists indeed, and the partial Hamiltonian at n = n c is 
defined on the domain of functions obeying condition 

lim( rMv/h) F f nctS (r) T(F) / @ x 

- * 2F ~ l 1 ] .tan(- + -), (50) 



\wi{rMv/h) l - F g ncS {r) Y{1 — F) \2 4 

r— >0 

where T(u) is the Euler gamma function, M is the parameter of the dimension of 
mass, and G is the self-adjoint extension parameter. Substituting the asymptotics 
of Eq.(49) at r — > into Eq.(50), one gets the relation fixing parameter u E , 

tm(u E ) = sgn(E) (^\ tan (| + |) ■ ( 51 ) 

In the case of graphitic nanocones with — | > rj > — 1 (jV^ = —4, —5) and 
f] < — 1 (Nj < —6), there are more than one irregular modes; this case will be 
considered elsewhere. 



4 Condensate 

It is instructive to rewrite Eq.(16) as 

p(x) = V-i(x), (52) 

where 

i(x) = -^^tr(x|aP|//| _1 |x). (53) 
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Although the trace of olP is formally zero, it may appear that current i is non- 
vanishing; then its nonconservation results in the emergence of condensate p. 

The contribution of regular modes is canceled upon summation over the sign 
of energy; thus, current (53) is vanishing in the case of 1 > 77 > |, and we are 
left with the cases of | > 77 > 0, > 77 > — |, and 77 = — 1 , when an irregular 
mode appears. Summing over s — ±1 corresponds to summing contributions of 
the inequivalent irreducible representations. These contributions are canceled for 
angular component ^(x) = —(i/A)hvtr{x\a v, P\H\^ 1 \x.) and doubled for radial 
component i r (x) = — (i/A)hv tr(x|a r P|i7| _1 |x). Consequently, we get 



i r U) 



x 




[L (+) + L H ] J^ F {kr) Ji- F {kr) + 



+ [L {+) - L H ][J F (kr) Ji- F {kr) - J_ F (kr)J- 1+F {kr)] 



hkV' 2F /e 7T 



cot hr + 7 + h-)]JF{kr)J- 1+F {kr) ) , (54) 



where 
L(±) = 



Mv V 2 4 



± tan — + - +2cos(Ftt)± — cot — + - 

[ Mv V 2 4 / y ' \Mv V 2 4 



-1 



(55) 

Extending the integrand in Eq.(54) to the complex fc-plane, using the Cauchy 
theorem to deform the contour of integration (for more details see Ref. [2]), and 
introducing the dimensionless integration variable, we recast Eq.(54) into the 
form 

00 

sin(FTr) r KfMKMw) f5 gx 

W ^(1-T]yj cosh[(2F-l)ln(^+lntan(f + f )] ' W 

where K^u) is the Macdonald function of order \x. Since in our case V • i = 
r~ x d r ri r , by differentiating Eq.(56) we get the following expression for the vacuum 
condensate: 

00 

sin(FTr) f K F (w) + K*_ F (w) 

Pl x ) = T - ^ / dWW p- ; 7 — e r 77^ rr. (57) 

K) n 3 (l-v)r 2 J cosh[(2F-l)ln(^)+lntan(f + f )] 1 ; 

Evidently, Eq.(57) vanishes at cosG = 0, while at F = | it is simplified, 

cos 

P X ) L_i = — — 7TT-, T^7- 58 

PV ;iF -2 2tt 2 (1 -r?)r 2 v ; 
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If cos 6 7^ and F ^ |, then at large distances from the defect we get 

. ,„ , ( frMv\2F-l r(f-F)r(f-2F) fQ ^ n^P^ 1 
sin(FTr) I {—) ffw) COt ( 2 + 4 J ' < b < 2' 



= _ Diil v^ J v n I t(i-f) 

- 2 (l-^ 2 \ m 1 - 2F ^to^tan(f + f), I<F<1. 



(59) 



5 Current 

It is straightforward to conclude that the radial component, j r (x) = 
— |tr(x|a r i?sgn(if)|x), is vanishing, so it remains to consider the angular 
component, J v (x) = — |tr(x|a ¥ 'i?sgn(iJ)|x). The contribution of irregular mode 
(49) to this quantity is 



V^'irregW 



47r(l - 1]) 



X 



X 



n v 



hk 



2F-1 



tan ( 7 + 7 ) [£(+)-£(-)U-F(fcr)Ji_ F (Ar) + 



+ [L (+) + L ( _)] [J F (£r) Ji_f(Ax) - J_ F (fcr) J_ 1+F (&r)] 

1-2F 



fife 



- ( 7 • T ) -i f _ ) ].7,,(^/-).7_ 1 _ / ,(A-/-) ^ . (60) 



where L(±) is given by Eq.(55). Similarly as in the previous section, we get 



7r 2 (l — ry)r s 



rfww< If(w)Ki-f(w) — Ii-f(w)Kf(w)+ 



2 sin(F7r) r ^ . w ^ . . 
+ -K F (w)K 1 _ F (w) tanh 



7T 



/ /ill) \ 



( 2F - 1 )" 1 (^J +lnta "i2 + 4 



71 



,(61) 



where I^u) is the modified Bessel function of order fi. The contribution of regular 
modes (47) and (48) is 



OO 

i f r 00 

n '=1 



yi ^'(l^)-i+F(^)^'(l-r,)-i-l+F(^) 



Z'=l 



(62) 
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Performing the summation (details will be published elsewhere), we get 



7r(l — r])r 2 



G(r],F) + - [F-- tan(F7r) + 



+ 



7T 



dww 



I F (w)K l _ F (w) - I l _ F (w)K F (w) 



(63) 



where 
G(r,,F) = 



= — I du- 

47T. 



sin(F7r) cosh 




+ sin 




cosh [(F-|) u] 


cosh 2 (f ) 


cosh (^) 


" cos (A) 





(64) 



Thus we get the following expression for the vacuum current: 



1 



7r(l — ??)r 2 

oo 

x J dw wK F (w)Ki- F (w) tanh 
o 



Gfa, F) + - F - - tan(FTr) - 



2 sin(FTr) 



x 



71 z 



(2F- l)ln 



hw \ 

rMv J 



. 7T 

+ In tan ( — + - 



(65) 



At cos = we get 



7r(l — ?])r 2 



G(t7, F) + -(1 ± 1) F - - tan(FTr) 



7T 



Also Eq.(65) at F = | is simplified, 



sin 



^ (x)| H = - 2^p - 



, = ±-(mod27r). 

(66) 

(67) 



If cos© 7^ and F 7^ |, then at large distances from the defect we get 



7r(l — r^)r 2 



C(77,F) + - F- 



) tan(F?r) 



(68) 



In the case of 1 > r\ > \ (N d = 5, 4, 3), the vacuum current takes form 



dk k 



1=1 



tt(1 -77) 

OO 

^ Ji'(i~ ri )-^+F{kr)J V (i- ri )-i-i +F (kr) 



^2 Ji(i- v )-i-F(kr)J l{1 _ ri) -i +1 _ F (k7 



(69) 
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In the cases of Nd = 3 and Nd = 4, the sums in Eq.(69) are canceled term by 
term; thus the current is vanishing. In the case of Nd = 5, the current can be 
presented in the following form 



Using the Schlafli contour integral representation for and K^iu), one can 

show (details will be published elsewhere) that the current is vanishing in this 
case also. 

6 Summary 

In the present paper we study the ground state polarization in graphene with a 
disclination, i.e. 6 — A^-gonal (Nd ^ 0) defect inserted in the otherwise perfect 
two dimensional honeycomb lattice. The variation of the bond length and the 
mixing of tt- with cr-orbitals caused by extrinsic curvature of the lattice surface 
are neglected, and our consideration, focusing on global aspects of coordination 
of carbon atoms, is based on the long-wavelength continuum model originating 
in the tight-binding approximation for the nearest neighbour interactions. Our 
general conclusion is that the ground state is polarized in cases when the Dirac- 
Weyl equation possesses a solution which is irregular, although square integrable, 
at the location of the defect; thus the ground state polarization is depending on 
the boundary parameter at this point, which exhibits itself as the self-adjoint 
extension parameter. The conclusion is consistent with the previously obtained 
result for the induced ground state charge in graphene with a disclination [3l H] . 

It is straightforward to demonstrate that the usual ground state current, 
(vacl^^a^lvac), and the ground state pseudospin-condensate, (vacl^E^lvac), 
are zero. In the present paper we consider other ground state characteristics: the 
P-condensate (9) and the i?-current (10), which in terms of the sublattice and 
valley field components (see Eq. (35)) are explicitly written as 




I 3 i +2 (w)K 3 i+i(w)] 



N d = 5. 



(70) 



p(x) 



(vac| *\ + (x)*b-(x) + W B +(x)* A -(x) + 
W A _(x)^ B+ (x) + *}j_(x)* A+ (x) \ |vac) 



(71) 



+ 



and 



y/93 v (x) 



(vac| -V a+ (x)*a-(x) + & b+ (x)*b-{x) 
V A _(x)V A+ (x) +¥ B _(x)V B+ (x)} I vac) 



(72) 
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(the radial current is vanishing). Whereas the current is invariant under time 
reversal, the condensate is invariant under time reversal and spatial inversion 
as well. In particular, in the chiral representation of the Dirac matrices (with 
diagonal 7 5 -matrix) one gets P = 7 and the condensate corresponds to the 
conventional chiral symmetry breaking condensate, (vac|\I/\I/|vac). 

In the cases of the one-pentagon (N d = 1), one-heptagon (N d = —1) and 
three-heptagon (N d = —3) defects, our results take form 



P(x) 



6 sin(7r/5) 
57T 3 r 2 



dww 



4/5 ^ 



™sh [|ln(^)-lntan(f + f )] ' 



2sin(7r/5) 

+ V 2 

N*=l, 



JL [ G _l tan(7r /5)_ 
57rr 2 | V 6' 5 / 20 V ' ' 



dw wK 1 / 5 (w)K i / 5 (w) tanh 



3 _ f hw \ . 

- in — mtan 

5 V rMv J 



N d = l, 

(73) 



e 71 
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P(x) 
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P(x) = - 
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dw w 
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— / dwwK 1/3 (w)K 2 /3(w) tanh 



iln ( , , 
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, 71 
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(78) 
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At large distances from the defect, the current decreases as r~ 2 , see Eq.(68), 
whereas the condensate decreases faster, see Eq.(59), with the same power law 
as for the decrease of the charge density [3]. 

In the cases of the two-pentagon (Na = 2), two- heptagon (Nd = —2) and six- 
heptagon (Na = —6) defects, the expressions for the condensate and the current 
are simplified and are given by Eqs.(58) and (67), respectively. Note that in these 
cases the charge is zero [3]. 

One can see that the ground state polarization effects cannot be eliminated 
at all by the choice of the value of the boundary parameter (O). Even in the case 
of cos© = 0, when the condensate and the charge are vanishing, the current is 
nonvanishing, see Eq.(66). The question of which of the values of G is realized in 
nature has to be answered by future experimental measurements, probably with 
the use of scanning tunnel and transmission electron microscopy. 
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